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BACKWATER EFFECTS OF OPEN CHANNEL CONSTRICTIONS 


H. J. Tracy’ and R. W. Carter,’ A.M. ASCE 


SYNOPSIS 


A method of computing the nominal backwater due to open channel con- 
strictions was sought. A practicable solution has been accomplished which is 
based on an empirical discharge coefficient and a laboratory investigation of 
the influence of channel roughness. Also investigated as a part of the labora- 
tory tests were effects of channel shape and constriction geometry. The so- 
lution invoives the computation of water-surface drop through the constriction 
and the determination of a factor which is the ratio of backwater to water 
surface drop. This ratio is shown to be a function of channel roughness, per- 
cent of channel contraction, and constriction geometry. 


INTRODUCTION 


It is generally recognized that the introduction of a constriction in an open 
channel causes local, and in some instances, far reaching changes in the water 
surface profile. Evaluation of this “backwater” effect is of interest to those 
concerned with the location of artificial obstructions in natural channels. 
Highway departments have found such information desirable as a part of hy- 
draulic studies made at proposed bridge sites. 

This study was undertaken as an extension of an investigation [1, 2]? re- 
cently completed by the U.S. Geological Survey, which was specifically con- 
cerned with the discharge characteristics of the constriction. Experimental 
data and certain computation procedures developed are directly applicable to 
the present study. Nomenclature common to the two investigations are used 
as defined for the earlier study. 


Object of the Investigation 

The object of this investigation was to study, by theoretical as well as by 
experimental means, the backwater effect of width constrictions in open chan- 
nels. The ultimate objective was a practicable method of computation adapted 
to the evaluation of this effect. 


Scope of the Investigation 

This study was limited to an investigation of single-opening constrictions, 
and to steady, tranquil flows. The hydraulic characteristics of the channel 
were assumed to be uniform over the backwater reach. Thus, in the labora- 
tory, the channels upstream and downstream from the constriction were made 
identical with respect to bed roughness, shape, and width. Variables in the 
laboratory investigation were discharge, depth of flow, channel roughness, 
and a limited variety of constriction geometries and channel shapes. It is 


1. Hydraulic Engineer, U.S. Geological Survey. 
2. Numbers in brackets refer to Appendix 11, References. 
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expected that the scope of this report will be extended to include multiple 
bridge openings as a part of continuing investigations. 
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Description of the Problem 


It is pertinent to a discussion of the backwater problem to compare the 
backwater condition with the normal, friction-controlled flow in the same 
channel reach without the constriction. 

The backwater effect due to a constriction in a channel of mild slope is 
first manifest upstream from the constriction by a backwater profile of the 
M-1 variety [3] in which the velocities and consequently the rate of loss of 
flow energy is less than the normal. The beginning of this reach is designated 
by section 0 in Fig. 1. Near the constriction, at section 1 in the same figure, 
the central body of water begins to accelerate, whereas deceleration occurs 
along the outer boundaries, and a separation zone forms in the corners up- 
stream from the constriction. 

At the constriction, the flow is characterized by rapid acceleration in di- 
rections both along and normal to the stream lines. The average longitudinal 
water surface profile falls rapidly in this region. Within the constriction, the 
live stream contracts to a width somewhat less than the gross width of the 
opening and the spaces between the live stream and the constriction boundaries 
are occupied by eddying water. The extent of the contraction and consequent 
separation depends upon rate of flow, constriction and boundary geometry, 
and roughness. 

Following the section of minimum live stream with (the vena contracta, 
section 2) is the beginning of an expansion process, which continues until the 
normal regime of flow has been re-established in the full width channel down- 
stream (section 4). The reach from section 2 to section 4 is one of decelerated 
flow in which the average velocities and energy losses are greater than nor- 
mal, due to the turbulence engendered by the expansion process. The length 
required to establish the normal regime downstream from the constriction is 
not known. Over the whole reach encompassed by the backwater effect, the 
total energy loss is the same as that for normal flow. 

As a specific objective, the theoretical and laboratory phases of this analy- 
sis aimed at a description of the backwater effect in a limited region of the 
“backwater reach.” Because of its use as a reference section in the earlier 
study [1], the nominal approach section (section 1, Fig. 1) was selected as the 
appropriate location for the measurement of the vertical extent of this effect. 

Section 1 has been defined as the beginning of acceleration of the flow ap- 
proaching the constriction. A distance upstream from the opening equal to 
one opening width (b) has been found to be an adequate approximation for the 
location of this section. Thus, the difference in water levels at section 1 be- 
tween the normal and backwater profiles is the backwater measure adopted. 
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This difference in water levels has been designated by the symbol hj on Fig. 1. 
Similarly, the backwater at the downstream sice of the constriction (section 3)5 


is denoted h3 on the same figure. 


* 


A convenient dimensionless measure of backwater is the ratio wr where 


Ah is the difference in water surface elevation between section 1 and the 
downstream side of the constriction, section 3. The definition of Ah in terms 
of the discharge equation and certain coefficients is taken from the earlier 
study [1], 


+ (1) 


where C, the coefficient of discharge, is a function of shape and roughness of 
the approach channel, an appropriate Froude number, and the geometry of 
the constriction. 

As illustrated in Fig. 1, 4h may be resolved into three components. From 
that figure, 


4b = & + 4h, 


where Ah, is the friction loss between sections 1 and 3 in the non-obstructed 
channel. This equation, rearranged to give an explicit expression for the 
nominal backwater, hj, becomes 


by = Ah bs - Ay 


and, in dimensionless form, 


1 


Ba ah dh 


Here, Ah is determined by the discharge characteristics of the upstream 
reach. For a given value of Ah, the ratio ‘3 is a measure of the relative 


magnitude of the recovery process, i.e., deceleration, and is a function of the 
downstream channel reach. The term - is representative of the part played 


by friction in producing Ah as compared with the influence of acceleration. 
The assumption is frequently made that the water surface immediately 
downstream from the constriction is the normal elevation, and that, neglect- 


ing friction losses, the backwater hf is equal to the drop, Ah. Accordingly, 


the backwater ratio be would then always have a value of 1.0. Equation (2) is 


useful in demonstrating that this assumption may be erroneous. The backwater 


3. Section 3 is arbitrarily defined. The piezometric level at section 3 is 
measured in the relatively quiet zones of eddying fluid at the downstream 
side of the constriction. The water surface level thus measured has been 
shown to be equivalent to the elevation of the water surface at section 2 [1]. 
The area at section 3 is computed using the average depth and gross width 
of the opening at the downstream side. 
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ratio is equal to unity only when the normal friction loss and the downstream 
recovery of energy are small in comparison with Ah. 


Approximate Solution for the Backwater Ratio 


An approximate analytical solution for the backwater ratio can be obtained 
if it is assumed that normal boundary friction losses are zero [4]. In the 
following development of the approximate solution it is also assumed that the 
channel is horizontal, that the velocities in the live stream at section 2 are 
essentially constant, and that the piezometric head at section 3 is equivalent 
to the level at section 2, as shown in the earlier study [1]. See Fig. 2. 

If Vo is the average velocity of the live stream at section 2, and V3 the 
average velocity at section 3 based on the gross area at section 3, the con- 
tinuity relationships are: 


Q = VyBy, = = V3by; = V,By, (3) 


From the momentum equation, 


wby3 w(B - 
2 2 


and from the energy equation, 


3g 


v2 
y: + + 
1 2g 2 


Using the equality yo = yg, and introducing the Froude number F3 = 
equations (3) and (4) may be combined to yield: 


3 23 & 2 
(2) i = 0 


% % 


Solved simultaneously, equations (3) and (5) yield: 


2 

F 
3 2+ b)2 


¥3 


Equations (6) and (7) for this approximate solution describe ¥3 and a 
¥4 3 


implicitly in terms of channel contraction ratio, 4 a Froude number, and 
the coefficient of contraction, C,. 

In the development [1] of an empirical expression for the coefficient of 
discharge C, that coefficient was shown to be primarily a coefficient of con- 
traction. A comparison of the expression for C with an equation derived theo- 
retically for the discharge through a constriction leads to the relationship: 
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2 
V3 


Ce 


VK 


in which K is used to describe eddy losses in the upstream separation zone, 
which are normally very small. The velocity distribution coefficient ao is 
commonly assumed to be unity. From dimensional analysis [1], C can be 
shown to depend upon a Froude number, the geometries of the upstream chan- 
nel and of the constriction, and upon an appropriate depth. 

Thus, C, may be represented by the coefficient C in equations (6) and (7). 
The variation of C with the variables named above is shown graphically in 
reference [2] in which C is defined for several practical constriction types. 

Two solutions for the backwater ratio are shown in Figs. 3 and 4. For 


these solutions, equations (6) and (7) were solved for ¥3 and 71 first by 


c = 


4 ¥3 
choosing constriction and channel geometries for which C = 1.00 for a full 
range of values of F and = A similar solution was made selecting C for 


constriction geometries having minimum contraction coefficients, those for 
vertical faced constrictions with square edged entrances. 


From these two solutions for 1 and ¥3, the backwater ratio for this ap- 
¥3 y 
proximate solution may be computed. For the frictionless condition, and with 
no bottom slope, att Hence: 


In Figs. 3 and 4 the description of the degree of channel contraction has 
been modified to be expressed as m = (1 - 7100 instead of the 2 ratio indi- 


cated in equations (6) and (7). 
For this special solution, it will be noted from figures 3 and 4 that the 
* 


Froude number has little effect on the ratio - and that this ratio is depend- 


ent primarily on m, with C, exerting a secondary influence. The effect of a 
* 


decrease in C, is to increase h1 but even when C, is varied through its 


Th’ 
maximum practical range, as it was in the two solutions, its effect is rather 
small. 


The Laboratory Investigation 


Exploratory tests in the initial phase of this study were made in the hy- 
draulics laboratory of the Georgia Institute of Technology. These tests es- 
tablished the relative importance of the variables involved and indicated the 
desirability of conducting further tests in a flume of greater length than that 
used for the exploratory tests. 
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Thus, arrangements were made for the use of a flume in the A. N. Talbot 
hydraulics laboratory of the University of Illinois at Urbana, Illinois. This 
flume is rectangular in cross section, 164 feet long, 5 feet wide, and 4.75 feet 
deep. The floor of the flume is level. A maximum discharge of 14 cfs is 
provided by the laboratory’s re-circulating systen:. Discharges are measured 
by means of a 12-inch bend meter in the supply pipe line. Baffles were used 
to obtain uniform velocities downstream from the supply entrance. 

The beginning of the test reach in the flume was about 15 feet from the 
supply pipe outlet. The depth of flow in the flume was regulated by means of 
an adjustable tailgate located 148 feet downstream from the beginning of the 
test section. Except in the vicinity of the constriction, water surface eleva- 
tions were read at 10 foot intervals along the length of the flume by means of 
fixed point gages set over the center of the flume. At the constriction, water- 
surface levels were read by means of a movable point gage which traversed 
the width as well as the length of the flume. Constrictive elements were of 
aluminum plate and plywood. For all tests, the constriction was located 
55 feet downstream from the beginning of the test section. 

Boundary roughness was created by means of various arrangements of 
chain-link fencing, and expanded metal. The cross-sectional shape of the 
test section was altered by means of false floors, using smooth transite sheets 
and redwood lumber. 


Methods of Testing 

The general test procedure for backwater measurements consisted of mak- 
ing two complete sets of surface-level observations over the full length of the 
test section. This was done, first, with the constriction in place and, second, 
with the constriction removed. 

Tests on the rectangular flume were made with four degrees of channel 
roughness. The concrete floor and sides of the basic flume provided a con- 
dition of minimum roughness corresponding to a Manning’s n of .012 which 
was constant over the range of depths tested. A second degree of roughness 
was created by covering the flume floor uniformly with 2-inch chain-link, 
diamond-mesh, wire fencing. The resulting n as determined in the laboratcry 
was 0.023. By adding expanded-metal sections which projected 2-1/2 inches 
above the chain-link fence, a value of n of .050 was measured at a flow depth 
of one foot. For this condition of roughness, n decreased slightly at greater 
depths, and was considerably greater at depths less than one foot. 

For testing the effect of a change in shape of channel cross section, a false 
bottom was built into the flume. The dimensions of this alternate arrange- 
ment are shown in Fig. 5. A limited number of tests were made in this flume, 
for which Manning’s n was computed as .010. 


Analysis of Experimental Data 


The laboratory phases of the investigation was first concerned with the 


measurement of the backwater ratio hy in a uniform rectangular flume at 


ah 
various degrees of bed roughness, using vertical-faced constrictions with 
square-edged abutments. Using the curves defined by this data as base, or 
standard values, other constriction and channel geometries were tested and 
the variation of the backwater ratio then determined as a departure from this 
base. 
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The Backwater Ratio for the Base Condition 
Sufficient tests were made to describe the backwater ratio for vertical- 

faced constrictions with square-edged abutments in the rectangular flume for 

three degrees of channel roughness, ranging from smooth to fairly rough. 

The results, plotted separately for each of the three roughness conditions, 

are shown in Figs. 6, 7, and 8. For convenience, these three curves have 

been summarized in Fig. 9, omitting the Froude number as a variable because 

of its minor effect. Also shown on this figure are curves resulting from the 

approximate analysis. It is apparent from the figure that the approximate 

solution is an aid in describing the shape and limits of the experimental curves. 
Important during this series were tests made to isolate the influence of 


depth of flow, 3 ; and of constriction length, é With all other variables held 


constant, a variation of these parameters between tests proved to have no 
effect on the backwater ratio. 


In this study, the width ratio 5 (or 1 - 4) has been sufficient up to this 


point to describe the percent of channel contraction. However, when irregular 
channels of non-uniform roughness are considered, this definition of the rela- 
tive degree of constriction imposed is no longer valid. By interpreting mas 
the proportion of total flow entering laterally into the contracted stream, a 
useful concept is gained. This may be stated thus: 


« $)100 


Where q is that portion of the total discharge Q occupying an area of width b 
in the total cross section upstream from the constriction. By computing the 


ratio 3 as a ratio of hydraulic conveyances, this definition may be stated: 


= (l- 2)100 
Q 


where Ky is the conveyance of that portion of the approach channel occupied 
by the discharge q, and is the conveyance of the total section. For the 
rectangular channel of uniform roughness characteristics heretofore con- 


sidered, this expression reduces to the previous of m, that is, (1 - =)100. 


In discussing equation (2), it was pointed out that the backwater ratio is 
influenced by the channel characteristics of both upstream and downstream 
reaches, as well as by the constriction geometry. It is reasonable to believe 
that the variable, m, descriptive of channel shape and roughness only at chan- 
nel section upstream from the constriction, is not necessarily always an ade- 
quate measure of channel geometry in the downstream reach. Actually, as 
tested in the laboratory, the channels were identical in shape and roughness 
in both reaches; therefore results developed here should not be applied to 
channels in which these two reaches are significantly different. 


The Effect of Channel Roughness 

Despite the limitations of Manning’s n as a measure of channel roughness, 
its use in this study is not a serious shortcoming, primarily because rough- 
ness is relatively unimportant as a factor in the determination of the backwater 
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ratio. This is shown in Fig. 10, which is an alternative presentation of the 
data shown in Fig. 9, replotted to show more clearly the variation of iz with 


n. This figure indicates that at an n of about .050 the limit of the change in 
the backwater ratio due to channel roughness has been reached. 

It is of interest to note that hj increases as n increases. This may be 
shown by an examination of equation (1) and Fig. 10. Equation (1) indicates 


an increase in Ah withn. Fig. 10 demonstrates an increase in pe with n. 


It is thus evident that hj must also increase with n. Similiarly, it may be 
shown that h3 decreases as n increases. 


The Effect of Channel Shape 

Twelve tests were made using the flume shown in Fig. 5 to determine the 
effect of a single variety of complex channel cross section on the backwater 
ratio. Vertical-faced constrictions with square-edged abutments were used 
as before to permit these tests to be compared directly with those made in 
the rectangular channel using the same type of constriction. Manning’s n 
for this channel was computed as .010. These twelve tests are plotted on the 
curves developed for the rectangular channel, and are shown in Fig. 11. Tak- 
ing into consideration the differences in channel roughness, the agreement 
between these tests and those made in the concrete rectangular flume (n = .012) 
is good. This has been taken as evidence that the influence of cross-sectional 
shape is incorporated in the variable m. 


The Effect of Constriction Geometry 

As suggested earlier, the effect of constriction type on the backwater ratio 
has been analysed as a systematic variation from the base curves of Fig. 9. 
Necessary to this treatment is an index descriptive of those features of geom- 
etry which influence flow phenomena in the backwater reach. It is reason- 
able to use the coefficient C for this purpose, since C is also primarily a 
function of constriction geometry. 

A number of tests were made in the rectangular, uniformly smooth con- 
crete flume using constrictions providing less abrupt entrances than those of 
the base type. For these tests, backwater ratios were obtained which were 
smaller than those shown by the base curves of Fig. 9. As a measure of the 
variation between the tests and the base curves, each of the backwater ratios 
of the tests was divided by the corresponding base ratio, selected from Fig. 9. 
The resulting factor, k,, is thus an expression for the departure of the tested 
condition from that of the base condition. 

By dividing the coefficient C applicable to the tested constriction by a co- 
efficient C appropriate to a vertical-faced, square-edged constriction with an 
identical m, a factor which correlates with k, is obtained. The variation of k, 
with this factor is shown in Fig. 12, in which the plotted points represent tests 
made on other than base type constrictions. In effect, this figure shows the 
variation in backwater ratio due to constriction geometry as a function of the 
coefficient of discharge. 


Computation Procedure 


This study has defined a nominal backwater from open channel constric- 
tions, which has been referred to as hj (refer to figure 1). This backwater is 
the difference between the normal and the constricted channel water surface 
profiles at an arbitrarily chosen section, section 1. This backwater index has 
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been divided by a known quantity, Ah, to give a dimensionless measure of the 
backwater effect. 

When integrated with an earlier study which investigated the discharge 
characteristics of open channel constrictions, and using the concepts and 
charts developed during that study, the results of the present analysis may 
be used to predict hj for constricted natural channels, within the scope of 
this study. This additional material is to be found in references [1, 2]. 


The backwater, hj, is computed as the product of at and h. The difference 


in piezometric heads between section 1 and 3 for the constricted channel, Ah, 
may be computed from equation (1). The solution of this equation depends 
upon a knowledge of the channel properties at section 1 and 3, and upon the 
evaluation of the coefficient C. ‘ 


The backwater ratio itself, ee has been shown by laboratory data to be a 


function primarily of degree of channel contraction. The influences of bed 
roughness and constriction geometry are secondary. Variables characteris- 
tic of the flow, as Froude number and depth, are largely unimportant in their 
effect on this ratio. | 


In application, = is selected from figure 9, appropriate to the bed rough- 
ness and degree of contraction of the channel. Compatible with the conditions 
under which the laboratory tests were made, this roughness must be repre- 


sentative of an average in the channel reaches both upstream and downstream 
from the constriction. Some fatitude is tolerable in the precision of the selec- 


tion, because of the relative insenstivity of mi with n. Computation of m for 


channel of non-uniform roughness is described in reference [2]. 


After z of the previous step has been chosen, it is adjusted for constric- 


tions geometry effect by the factor k,, obtained from figure 12. This correc- 
tion factor is a function of a & Leanne the computation of which has been out- 
lined earlier, and the degree of contraction m. 

The adjusted backwater ratio may then be multiplied by Ah to give hj. 


APPENDIX I 
NOTATION 


area of channel cross section 
width of opening 

width of channel 

coefficient of discharge 


coefficient of contraction; ratio of area of contracted stream to area of 
opening 


Froude number; a dimensionless flow parameter for gravity-controlled 
phenomena 
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3. 


acceleration due to gravity 
piezometric head; sum of elevation and pressure heads 


a difference in piezometric heads involving normal and constricted 
channel flow 


conveyance; the quantity is ARS in Manning’s formula for open chan- 
nel flow. Also, an energy loss coefficient 


an adjustment coefficient; ratio of backwater ratio to the backwater 
ratio for the base condition 


length of constriction, measured parallel to flow 
channel contraction ratio 


roughness coefficient in the Manning formula for open channel flow; 
also used as a subscript to denote “normal” 


volume rate of flow; total discharge at a channel cross section 

average velocity; ratio of total discharge to total area of a cross section 
specific weight; weight per unit volume 

depth of flow referred to average bottom elevation of cross section 

a symbol denoting “change” or “difference” 


density; mass per unit volume 
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FIG.1- DEFINITION SKETCH OF THE BACKWATER REACH 
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FIG. 5- CROSS SECTIONAL VIEW QF FLUME ARRANGEMENT (ILLINOIS) 
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FIG.6- MEASURED VALUES. OF THE BACKWATER RATIO, = 0.012 
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FIG. 7— MEASURED VALUES OF THE BACKWATER RATIO, n= OC23 
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FIG.8- MEASURED VALUES OF THE BACKWATER RATIO, n= 0.050 
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FIG. 9~ SUMMARY CURVES OF THE BACKWATER RATIO FOR VERTICAL FACE CONSTRICTIONS WITH SQUARE EDGED ABUTMENTS 
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FIG 1O- THE EFFECT OF CHANNEL ROUGHNESS ON THE BACKWATER RATIO, VERTICAL 
FACE CONSTRICTIONS WITH SQUARE EDGED ABUTMENTS 
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FIG. \1- THE EFFECT OF CHANGE IN CHANNEL CROSS SECTION ON THE BACKWATER RATIO, VERTICAL FACE 
CONSTRICTIONS WITH SQUARE EDGED ABUTMENTS 
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FIG. 12— THE EFFECT OF CONSTRICTION GEOMETRY ON THE BACKWATER RATIO 
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